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The remarkable discrepancy between theory and data on the pt distribution of J/ip polarization 
in the process pp J/tp + X has been standing as a puzzle for a long time. Our early work 
indicates that if only considering the color-singlet state, NLO corrections may drastically change 
the distribution, but are not sufhcient to accommodate the data. The color-octet intermediate 
state was suggested to make up the gap. Therefore a direct motivation is to evaluate if the NLO 
corrections to the J/tp production via a color-octet can indeed remedy the problem. Unfortunately, 
our results show that even the NLO corrections to both color-singlet and octet states are accounted, 
the data cannot be understood yet in the framework of the present theory, but there, as observed, 
still a small window remains, namely a P-wave octet might substantially contribute. 
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Since its discovery in 1974, study on J/tp production 
never ends, however it seems that this field is still not 
fully understood. To solve the large discrepancy be- 
tween the experimental data and theoretical predictions 
on the Pt distribution of J/tp production at Tevatron, a 
color-octet mechanism [ij which may increase the large 
Pt distribution was proposed based on the non-relativistic 
QCD (NRQCD) [2]. However, despite the developments 
and successes, the predictions of NRQCD do not al- 
ways make satisfactory predictions for some processes. 
The authors of Ref. [3] find that the DELPHI '4\ data 
for J/tJj production in 77 — > J/^X evidently favor the 
NRQCD formalism, but a certain amount of color-octet 
must be included. Furthermore, an analysis on char- 
monium production at fixed-target experiments was per- 
formed within NRQCD up to NLO [8j] , and it is indicated 
that fraction of color-octet which is needed to describe 
the data is only about 1/10 of that required for explaining 
the Tevatron experiment. By contrast, the experimental 
results on inelastic J/4> photoproduction at the DESY 
ep collider HERA are adequately described by the color- 
singlet mechanism alone once higher-order QCD correc- 
tions are included I,!]. At B-factories, JMfPlf^] pro- 
duction suggested by the authors of Ref. was not dis- 
covered. Even more seriously, the LO NRQCD calcu- 
lation predicts a sizable transverse polarization rate for 
large pt J/jt whereas the Tevatron measurement at 
Fermilab [l3| displays a slight longitudinal polarization 
at large pt ■ 

On the other hand, obvious discrepancies between LO 
predictions [ll|, and experimental results [H, [l3| for 
single and double charmonia productions at B-factories 
are remarked. Further studies indicate that they may be 
resolved by including higher order corrections: both NLO 
QCD and relativistic corrections [UlIlBIlBl, or at least 
alleviated. Recently, the NLO QCD corrections to J/ip 
hadronproduction are calculated in Refs. [17,, ilSj] . The re- 
sults show that the production rate of J/tp at the larger 
transverse momentum pt region is much increased. The 



NLO process gg — > J/tpcc is calculated in Refs. [H, 
and it causes a sizable contribution to the pt distribu- 
tion. The NLO QCD corrections to J/ip polarization via 
color singlet at Tevatron and LHC have been calculated 
in Ref. [18| and the results show that the J/ip polarization 
status drastically changes from transverse-polarization 
dominance at LO into longitudinal-polarization domi- 
nance at NLO, however, such changes still cannot ex- 
plain the data yet. As the NLO corrections are so im- 
portant for the Pt distribution and polarization status 
and as supposed, the color-octet mechanism is needed 
to explain data, it is obviously necessary to investigate 
whether NLO corrections can seriously affect J/ip pro- 
duction via color octet states. In this letter, we calculate 
the NLO QCD corrections to the J/tp production via 

color octet states J/tp[S^\ at Tevatron and LHC 

where the Feynman Diagram Calculation package (FDC) 

[2l| is employed. However, it is noted that J/^['^j*-'] is 
not included since the part for dealing with P-wave loop 
processes in FDC is not completed yet. 

According to the NRQCD factorization formalism, 
the inclusive cross section for direct J/ip production in 
hadron-hadron collision is expressed as 

a[pp^J/ip + X] = J2 J dxidx2G,/pGj/p 

xa[z-f j^(cc)„-fX](0,f), (1) 

where p is either a proton or an antiproton, the indices 
i,j run over all the partonic species and n denotes the 
color, spin and angular momentum states of the inter- 
mediate cc pair. The short-distance contribution a can 
be perturbatively calculated order by order in as- The 
hadronic matrix elements (O^) are related to the tran- 
sition probabilities from the state (cc)„ into J/ip which 
are fully governed by the non-perturbative QCD effects. 
In the following, a represents the corresponding partonic 
cross section. 

At LO, there are three partonic processes: 
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J/^[^'\¥,%ps) + q{p,h (L2) 
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g{Pi) + q{P2) 
q{pi) + q{p2) - J/i^['S^'\"sf%3)+9{P4). (L3) 



where q represents a sum over all possible light quarks or 
anti-quarks: u, d, s,u, d, s. 

The NLO corrections include virtual and real cor- 
rections. There exist UV, IR and Coulomb singulari- 
ties in the calculation of the virtual corrections. UV- 
divergences from self-energy and triangle diagrams are 
removed by renormalization. Here we adopt the renor- 
malization scheme used in Ref. [22j |. The renormaliza- 
tion constants Z2, Z21 and Z3 which correspond to 
charm quark mass mc, charm- field light quark field ipq 
and gluon field A'^ are defined in the on-mass-shell(OS) 
scheme while Zg for the QCD gauge coupling as is de- 
fined in the modified-minimal-subtraction(MS) scheme: 
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where 7f is the Euler constant, /3o = — ^TpUf is 

the one-loop coefficient of the QCD beta function and 
Uf = 3 is the number of active quark flavors, is the 
renormalization scale. 

There are 267 (the ^S^^^ state) and 413 (the ^ state) 
NLO diagrams for process (|L1[) . including counter-term 
diagrams, while for both processes (|L2[) and (jL3p . there 
are 49 (for the state) and 111 (for the ^S[^^ state) 
NLO diagrams altogether. Diagrams where a virtual 
gluon line connects the quark pair possess a Coulomb 
singularity, which can be isolated and attributed into the 
renormalization of the cc wave function. 

For each process, by summing over contributions from 
all diagrams, the virtual corrections to the differential 
cross section can be expressed as 



dt 



cx 2Re(Mf M^*) 



(3) 



where Mf is the amplitude of process (i) at LO, and 
is the renormalized amplitude of corresponding process 
at NLO. is UV and Coulomb finite, but still contains 
IR divergences. 

One also notices that as the real corrections where a 
real gluon is emitted from either initial or final quarks 



are taken into account, the three subprocesses (LI), (L2) 
and (L3) tangle together and must be considered simul- 
taneously. There are eight processes involved in the real 
corrections: 
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where q' denote light quarks (anti-quarks) with dif- 
ferent fiavors. Phase space integrations of above pro- 
cesses generate IR singularities, which are either soft or 
collinear and can be conveniently isolated by slicing the 
phase space into different regions. Here we adopt the 
two-cutoff phase space slicing method [2^ to deal with 
the problem. Then the real cross section can be written 
as 
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It is observed that the IR singularities from one real pro- 
cess may be factorized into different parts and each of 
them should be added into the cross sections of different 
LO processes. This is the reason why we have to calculate 
the NLO corrections to the three LO processes together. 

from the soft region contains soft singularities and 
is calculated analytically under the soft approximation. 
One should notice that, unlike color singlet case, the soft 
singularities caused by emitting a soft gluon from the 
charm quark pair in the S-wave color octet exist and 
the factorized matrix element is the same as the case 
where a soft gluon is emitted from a gluon. a^'~^ from 
the hard collinear region contains collinear singularities 
which are factorized out and the singularities are partly 
absorbed into redefinition of the parton distribution func- 
tion (PDF) (usually called as mass factorization [24i|). 
Here we adopt the scale dependent PDF using the MS 
convention given by Ref. [2J]. After redefining the PDF, 
an additional finite term <J^^ is separated out. The hard 

noncoUinear part a"^ is IR finite. Finally, all the IR sin- 
gularities are canceled and -\- a^'^ + is IR finite. 

To obtain the transverse momentum pt distribution of 
J/V', a transformation of integration variables {dx2At — > 
JdptAy) is needed. Then we have 
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where y is the rapidity of J/?/' in the laboratory frame and 
/i J is the factorization scale. The polarization parameter 
a is defined as: 



a{pt) 



dar/dpt - 2daL/dpt 
dar/dpt + 2daL/dpt 
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To evaluate a{pt), the polarization of J/iJj must be explic- 
itly retained in the calculation. The partonic differential 



FIG. 1: Total cross section of J/ip hadronproduction at LHC 
(upper curves) and Tevatron (lower curves), as function of /i 
with fir = fJ'f ~ l-t. and = \/ (2mc)^ + p^. The center-of- 
mass energy are 1.98 TeV at Tevatron and 14 TeV at LHC. 
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FIG. 2; Transverse momentum distribution of J/^j) produc- 
tion with fir = fJ-f = l-io at LHC (upper curves) and Tevatron 
(lower curves). 

cross section for a polarized J/i/j is expressed as: 
^ = ae{\).e*{X)+ ^ a,, ■ eiX) p, ■ e* (X), (7) 

where A = ri,r2,L. e{Ti), e{T2), e(L) are the two trans- 
verse and longitudinal polarization vectors of J ftp respec- 
tively, and the polarizations of all the other particles are 
summed over in n-dimension. One can find that a and 
CLij are finite when the virtual corrections and real cor- 
rections are properly handled as aforementioned. The 
gauge invariance is explicitly checked by replacing the 
gluon polarization vector into its 4-momcntum in the fi- 
nal numerical calculation. 

In our numerical computations, the CTEQ6L1 and 
CTEQ6M PDFs ^2^ are used for LO and NLO calcu- 
lations respectively. The charm quark mass is set as 
1.5 GeV. By fitting the pt distribution of prompt J/ip 
production measured at Tevatron [2^, the NRQCD 

matrix elements (O^) are determined as (Og (%*!)) = 
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FIG. 3: Transverse momentum distribution of prompt J/tp 
production at Tevatron. CDF data is from ref |26j]. 



0.0045 GeV^ and (OtCSo)) = 0.076 GeV^ and the re- 
sults are shown in Fig. [3l In the fitting procedure, con- 
tributions from both color singlet and octet states are 
included. However, it is worth noticing that we have 
to abandon the experimental data with pt < 6 GeV, 
since it is impossible to obtain a satisfactory pt distri- 
bution in terms of a unique (O^) value. In addition, 
one should consider an additional contribution of the 
feed-down from -0' which may bring up an extra factor 
B{^P' J/ijj + X) X {Oii)/{Ot), a short calculation de- 
termines it as 1.29 for color singlet part, and the values 
of our fitted {OtCSi)) and (O^(^5o)) include the contri- 
butions of the feed-down for color octet part. However, 

the feed-down from and Xcj is not considered 

in our fitting since at NLO it cannot be properly calcu- 
lated so far and this omission should be treated as an 
approximation. The two phase space cutoffs 5s = 
and 5c = (5^/50 are chosen, and the invariance for differ- 
ent values of 5s and 5c is obviously observed within the 
error tolerance. The result is restricted to the NRQCD 
applicable domain > 3 GeV, and \yj/^\ < 3 for LHC, 
\yj/ip\ < 0-6 for Tevatron respectively. 

The dependence of the total cross section on the renor- 
malization scale fir and factorization scale fif are shown 
in Fig. [TJ It is obvious that the NLO QCD corrections 
make such dependence milder. The pt distributions of 
J/^ production are presented in Fig. [2] where only slight 
change appears when the NLO QCD corrections is in- 
cluded. J/ip\^S(p] produces unpolarized J/ip for both 
LO and NLO. The pt distributions of J/4> polarization 
parameter a for J/^/; ['^S'^^^] are shown in Fig. 3] and there 
is a slight change when the NLO corrections are taken 
into account. 

As a summary, in this work, we have calculated the 
NLO QCD corrections to J/4> production via color octet 
states J/ij[^S^^\^s[^^] at Tevatron and LHC. With fir = 
fif = /io , the K factors of total cross section (ratio of NLO 

to LO) are 1.235 and 1.139 for J/i}[^S^^^] and J/i/'[%'f ^] 
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FIG. 4: Transverse momentum distribution of polarization a 
for prompt J/i/> production, CDF data is from ref [lol |. 

at Tevatron, while at LHC they are 0.826 and 0.800 re- 
spectively. Unlike for the color singlet case, there are 
only slight changes to the transverse momentum distri- 
butions of J/ip production rate and the J/ip polarization 
when the NLO QCD corrections are taken into account. 
The results imply that the perturbative QCD expansion 
quickly converges for J/^p production via the S-wave color 
octet state, in contrast with that via color singlet, where 



the NLO contributions are too large to hint a convergence 
at the NNLO. Finally, a definite conclusion is made that 
comparing the theoretical results calculated up to NLO 
with the experimental measurements at Tevatron, obvi- 
ous gap is observed, even though both color singlet and 
octet are included. In the well established theoretical 
framework of NRQCD there still remains a narrow win- 
dow which might make up the gap, namely one needs to 
investigate the NLO corrections to J/tp production via 
P-wave color octet state and J/ijj production by feed- 
down from Xcj- It is unclear how the situation will be 
when contributions from these two sources at NLO are 
taken into account, as we know that NLO QCD correc- 
tions to P-wave state in e+ + e^ J/i'+XcO evaluated in 
Ref. [27*1 are very large. However, if we assume that the 
NLO QCD corrections to the two P-wave states are just 
as those for the S-wave color octet states, it is reasonable 
for us to reach the conclusion that the large discrepancy 
of J/ip polarization between theoretical predication and 
the experimental measurement cannot be solved by just 
including NLO QCD corrections within NRQCD frame- 
work and new solutions should be explored. 
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